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Information about the pairing mechanism for superconductivity is contained in the spectral weight 
for the anomalous (Gorkov) Green function. In the most general case, this spectral weight can change 
sign on the positive real axis or even be complex in the presence of broken time-reversal symme¬ 
try. This renders impossible the direct analytic continuation with Maximum Entropy methods of 
numerical results obtained for the anomalous Green function either in Matsubara frequency or in 
imaginary time. Here we show that one can define auxiliary Green functions that allow one to 
extract the anomalous spectral weight in the most general multi-orbital spin-singlet or triplet cases 
with no particular symmetry. As a simple example, we treat the case of lead in Eliashberg theory: 
Maximum Entropy analytic continuation with auxiliary Green functions (MaxEntAux) allows us to 
recover spectral functions that change sign on the positive real axis. The approach can be extended 
to transport quantities. 

PACS numbers: 74.20.-z, 74.20.Mn 


Not so long after the proposal of the phonon mech¬ 
anism for superconductivity in BCS theory [1, 2], the 
central role played by retardation was recognized [3, 4], 
The convincing proofs of the validity of the phonon pair¬ 
ing mechanism, implicit in BCS theory [5], came from 
small deviations from BCS explained by the more re¬ 
fined theory of Migdal and Eliashberg [4, 6] . One of the 
most direct proofs came through the frequency depen¬ 
dence of the gap function measured through tunneling 
experiments [7]. Indeed, signatures of phonon frequen¬ 
cies and their harmonics observed in the frequency de¬ 
pendence of the density of states were explained in detail 
by theory [8-10]. 

In the search for pairing mechanisms in unconven¬ 
tional superconductors, such as cuprates and organic su¬ 
perconductors, one must face the presence of strongly- 
correlated physics, manifest through the proximity be¬ 
tween the Mott insulating and superconducting phases. 
The discussion between P. W. Anderson [11], and D. J. 
Scalapino [12, 13] shows that strongly-correlated physics 
sets again the question of retardation, or more generally 
the frequency dependence of pairing, at the heart of the 
problem. New high-energy pairing mechanisms could be 
involved. 

The thorough theoretical study of this subject is not 
an easy one as it requires to have access to the frequency 
dependent pairing dynamics embodied in the anomalous 
self-energy or in the Gorkov function [14]. While direct 
real-frequency studies are sometimes possible [15-18], 
most calculations, especially at finite temperature, rely 
on Matsubara-frequency or imaginary-time calculations. 
Direct analytic continuation with Pade approximants is 
occasionally possible [19, 20], but for Quantum Monte 
Carlo data in particular [21, 22], one must rely on Max¬ 
imum Entropy methodology for analytic continuation to 


the real axis [23]. However, this method works only for 
spectral functions whose sign does not change over the 
positive or negative frequency range. The anomalous self¬ 
energy and spectral function for the Gorkov function do 
not satisfy this criterion. A workaround using a constant 
shift has been suggested [24], but this involves another 
adjustable parameter and is not valid for complex spec¬ 
tral functions. 

We first consider the simple case of singlet pairing in a 
one-band model and give an example of application with 
the gap function for lead where the answer is known. The 
more general case is considered at the end. 

Auxiliary Green Function for Analytic Continuation 
of T Green functions associating a creation or annihi¬ 
lation operator to its hermitian conjugate do have a con¬ 
stant sign spectral weight and can be analytically contin¬ 
ued using Maximum Entropy methodology. With this in 
mind, following Ref. [21], define the mixed operator 


and the auxiliary Green function 

g^(k,r) = -(f T a^r)a * (0))^. ( 2 ) 

where the average is with respect to the grand canoni¬ 
cal ensemble for Hamiltonian TL and the imaginary-time 
ordering operator is 1~ T . In terms of the original annihi¬ 
lation (creation) operators c^, the auxiliary Green func¬ 
tion is given by 

S“(t,T) = -{ti i ,Wit t (0)) i - (teV (T)C^fO))^ 

- ('tej t W£_ J j(0)) s - ('te , _ ft H£i.(0)} s (3) 
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If the maximum-entropy analytic continuation is 
done from imaginary time, it suffices to analyt¬ 
ically continue both Q aux (k,r) and Q^{±.k, r) = 

— ( T t c , (t)c , (0)) that have positive spectral 

\ ±fct4- / it 

weight and to extract the spectral weight for the off- 
diagonal terms (Gorkov functions) by taking the differ¬ 
ence. This is possible because in the case considered, 
the Gorkov function F(k,r) = — (t t c^(t) c_^(0)^) _ 

and its complex conjugate have identical spectral weight. 
Indeed, for singlet pairing with time-reversal symmetry, 
the spectral weight A an (k,oj) is real and odd in fre¬ 
quency, implying that the Gorkov function in Matsubara 
frequency is also real 


F{k,iw n ) = / 7 T“ 


d oj A an (k,u) 

2n 


lUJ n — U! 

'dw a iA an ik,uj) 
27r u ; 2 + w 2 


( 4 ) 


Hence, the Matsubara frequency result 


Q aux {k, iu> n ) = G' t {k,iu n )-Gi(k,-iuj n )+2J r (k,iu n ) (5) 


immediately translates into 


A an {k,oj) 


A aux (k,oj) - A^(k,oj) - A±(k, -w) 


(6) 

with no restriction to the particle-hole symmetric case, as 
assumed in Ref. [21]. The spectral weight for the auxiliary 
Green function A aux (k,u}) is normalized to two instead 
of unity. 

An Example with Sign Changing Gap Function Let us 
focus on the analytic continuation of the solutions to the 
Matsubara frequency Eliashberg equations for lead. Usu¬ 
ally, this analytic continuation is performed with Pade 
approximants. Here, we take published real-frequency 
results [10] that we transform to Matsubara frequency. 
Our method, MaxEntAux, is then tested by analytically 
continuing these to real frequencies and comparing with 
the original data. 

Consider the spline interpolation A(w) of the gap func¬ 
tion extracted from Ref.[10], shown in Fig.l. The trans¬ 
verse phonon frequency w* = 4.4 meV is taken as unity. 
From this gap function, one can work out the local and 
anomalous spectral functions 


A 1oc (oj ) =CRe 


A an {oj) = C Re 


(v^ 2 -A 2 (u;)) 

( ah ^ 

\ v^ 2 ~ A 2 (w) y 


( 7 ) 

( 8 ) 


where the square root is chosen such that the real part is 
odd in oj and the imaginary part even, C is the constant 
normalizing A 1oc (uj) to unity over the frequency interval 
[— u! c ,w c ], where the cutoff frequency is to c = 25. Note 


that since the published real and imaginary parts of A(w) 
extend only up to w ~ 7, the resulting spectral functions 
(7) and (8) are extended as constants up to uj c such that 
both real and imaginary parts of A(w) obtained from 
Kramers-Kronig relations and from Eq.(13) agree with 
the published results [10]. 
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Figure 1. Solid magenta and green lines : real and imaginary 
parts of the gap function extracted from Ref. [10]. Black solid 
lines : spline interpolations of the data. Note that these inter¬ 
polations are extrapolated within the gap with a polynomial 
function that satisfies the continuity of the gap function and 
its derivative, and has a zero slope at u> = 0. Im A is odd in 
uj and Re A even. 


Using the spectral representations 



dw A loc (u) 

2i r iu> n — w 


F(iw n ) 


'dw A an (uj) 


2n icon — w ’ 


(9) 

( 10 ) 


Eqs. (7) and (8) give the local and anomalous Green func¬ 
tions in Matsubara frequency for T = 1.3 K = 0.02 w\. 
The integral over wave vectors of Eq.(5) yields 


G aux {iu n ) = G loc {iu n ) - g loc (-ioj n ) + 2F{ico n ). (11) 


We add a gaussian noise to impose a relative error 
of 10~ 4 on each Matsubara frequency component of the 
functions in Eq.(ll) to simulate Monte Carlo statistical 
errors. [25] Using the Maximum Entropy code Omega- 
MaxEnt [26], we then analytically continue g loc (iL 0 n ) 
and g aux {iu> n ) to real frequency to obtain A 1oc (oj) and 
A aux ( oj). Any reliable MaxEnt code should give similar 
results. The anomalous spectral weight is then extracted 
from the integral over wave vectors of Eq. (6) 

A an (co) = i [A aux (u>) - A 1oc {uj) - A loc (- w)] . (12) 

The imaginary parts of the local and anomalous Green 
functions are proportional to A loc (u}) and A an (uj) while 
their real parts are obtained using Kramers-Kronig rela¬ 
tions. The gap function then follows from [27] 


A(w) 


2wT(u) 

G loc {u) - g loc *(-uj ) ‘ 


(13) 
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Gap function One can now compare the original gap 
function of Fig.l with the gap function Eq.(13) obtained 
from maximum-entropy analytic continuation. A direct 
comparison is shown in Fig. 2. 



Figure 2. Red solid lines : spline interpolations of the real and 
imaginary parts of the gap function extracted from Ref. [10], 
before performing the maximum-entropy analytic continua¬ 
tion. Blue solid lines : real and imaginary parts of the 
gap function obtained from Eq.(13), after performing the 
maximum-entropy analytic continuation. 


herence length, the frequency dependence at the same 
position informs us on the dynamics. 




Figure 3. Red solid lines : anomalous spectral function 
and cumulative order parameter extracted from Eq.(8), be¬ 
fore performing the maximum-entropy analytic continuation. 
Blue solid lines : anomalous spectral function and cumula¬ 
tive order parameter obtained from Eq.(12), after performing 
the maximum-entropy analytic continuation. Dashed lines : 
frequency of the maximum of each cumulative order param¬ 
eter. Green dotted line : value of the superconducting order 
parameter. 


Even though some of the details are lost in the analytic 
continuation, the main features are preserved, namely 
the gap function has a maximum and becomes negative 
at sufficiently high frequency before vanishing. 

Gorkov function and Cumulative order parameter 
The top panel of Fig. 3 compares the anomalous spectral 
function obtained from Eq.( 8 ) using the published A(w), 
with the same function extracted from Eq.(12) after an¬ 
alytic continuation of the Matsubara data. A negative 
component is clearly visible in both quantities above 2.5 
times the transverse phonon frequency. 

The cumulative order parameter, exhibited in the bot¬ 
tom panel of Fig. 3, contains the important information 
about the dynamics of pairing and it is clearly well pre¬ 
served by the analytic continuation. At finite tempera¬ 
ture, it is defined by: 

iAu}) = (i4) 

where f(-uj') = [1 + e~^ u ] _1 is the Fermi-Dirac distri¬ 
bution. In the limit u> —> oo, we find that this quantity 
is just the order parameter Ijf(oo) = (%■ c if cq,^ _. At 
T = 0, Eq.(14) reduces to the formulas in Refs. [16, 17]. 
In the same way that the dependence of (%- Cjij . 
on distance i — j tells us about the superconducting co- 


As can be seen form Fig. 3, the frequencies larger than 
the gap contribute positively to the order parameter 
whose asymptotic value is represented by the horizontal 
dotted line. All frequencies up to the maximum, around 
2.5 times the transverse phonon frequency, contribute to 
pairing. After that, the effect of the Coulomb repulsion 
represented by fi* in the Eliashberg equation are depair¬ 
ing. The position of the peak (vertical dashed line) and 
the region where depairing occurs are well preserved by 
the analytic continuation. 

Auxiliary Green function for Analytic Continuation for 


General Multi-band Singlet or Triplet Pairing 

Define, 

Gj(k, t) = - (Tt%(t)c^(0)^ 

(15) 

where 7 stands for both band and spin index 
define the mixed operator 

. Again, 

^’k'yS ~ kS ’ 

(16) 

and the auxiliary Green function 



(17) 


The Lehmann representation allows one to check that 
this Green function has a positive spectral weight and can 
thus be analytically continued using standard Maximum 
Entropy methods. 

From the imaginary-time definition, one finds 
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OT'&t) = ^(k,r) - G«{-k, -t) (f T c^(r)c_^(0)) n - (f T c^Jr) 6^(0))^ . (18) 

The Lehmann representation of each of the last two terms, 

- (^%(t)c_£ 4 (0))^ = -^^e- pH ™{m\e H ™ T c h e- H ™ T \n){n\c_ %s \m) (19) 

mn 

- (tr c f _^(r) 4 7 (°))^ = E e ^ m He^c^e-^InXnlc^lm), (20) 

mn 

where Z is the grand partition function, and the equality [{m\ct \n){n\ct s \rn)]* = (m\c^ - |n)(n|ct |m) show that in 
general the spectral weights of the two anomalous functions are complex conjugate of each other. Using the spectral 
representation of the anomalous spectral function, 


F l5 {k,iuj n ) = ~ Are 1 ' 


% 7 ( r ) c_g 5 ( 0 )— / 0 _ 


'dw A$(k,w) 


n 


27t iw„ — w ’ 


Eq.(18) can then be written as follows in Matsubara frequencies 


g*rHk,iun) = Gj(k,ioj n ) - Gs{-k,-iu n ) + 


Aco A^(k,u) + A^*{k,u) 

2n iui n — ui 


( 21 ) 


( 22 ) 


In practice, the auxiliary Green function in Matsubara frequency, that we need to analytically continue, would be 
computed using the results of Matsubara frequency calculations as follows: 

G^ xl (k,ioj n ) = G y (k,iuJn) - Gs(~k, -iw n ) + F^s(k,iUn) + [^s(k, -ioj n )}*. (23) 

Eq.(18) can also be used instead if the analytic continuation is performed directly from imaginary time. 


To find the missing information, we define, by analogy 
with Eq.(16), a mixed operator 


and a corresponding second auxiliary Green function 
with positive spectral weight 



'■'k-y 


i c 


- kS ’ 


= - (^w^U 0 ))* ■ (25) 

(24) Then, by analogy with Eq.(22), we find, 


G“s x2 (k,iu n ) = G~/(k, iu> n ) - Gs{~k , -iui n ) - i 


do; A^(k,uj)-A;r(k,u) 
27r ioj n — u 


(26) 


The analytic continuation of the two auxiliary Green functions gives the corresponding real and positive spectral 
weights, which are formally given by 

A a 7 T 1 (k,uj) = A^k,uj) + As(-k, -w) + AS$ (fc, w) + A%*(k, w) (27) 

A;r 2 (k,u) =A 1 (k,u) + A s (-k,-u J )-i(A^(k,u J )-A;r(k,u)). (28) 

The normal spectral functions are also real and positive. So, from the above, one can easily extract the needed spectral 
weight A“g(k,w) via 


A$(k,u) = \ [a“ u s x \ k,u) +iA“T 2 (k,u) - (1 + *) (Ay(k,u) + A s (-k,-bo)) 


(29) 


Recalling that 7 and S contain both spin and band index, 
this is valid for singlet or triplet pairing and for the multi¬ 
orbital case, with or without inversion or time-reversal 


symmetry. 

Conclusion The above results open the way to the 
systematic exploration of the pairing mechanism in nu- 
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merical calculations for unconventional superconductors. 
Similar auxiliary functions can be defined to find, with 
Maximum Entropy methods, spectral functions of trans¬ 
port coefficients, such as the thermoelectric power, that 
do not have positive spectral weight. 
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